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Abstract. We derive an explicit method for reconstructing singularities of the 
initial data in a thermoacoustic tomography problem, in the case of variable sound 
speed and limited boundary data. In order to obtain this explicit formula we assume 
the metric induced by the sound speed does not have conjugate points inside the 
region to be observed. 

1. Reconstruction of the Wavefront Set Given Measurements on a 

Hyperplane 

Let fi be a bounded region in R™ such that x n < for x G ft. Let H be the 
hyperplane {x n = 0}. We study the Cauchy problem 

Pu = u u — c 2 (x)Au = 

u(x,0) = f(x) 

(1) u t (x,0) = 0, 

where c is a smooth function which satisfies -h < c(x) < M for some number M > 
1. This is the Cauchy problem commonly studied in mathematical investigations of 
thermoacoustic tomography and photoacoustic tomography (see KuKy for a survey 
of results in this area). 

The differential operator P has symbol 

(2) p(x,tr)=r 2 -c 2 (x)\e . 

c(x) induces a Riemannian metric on R™ by setting gij(x) = c~ 2 (x)5ij. We assume 
c(x) equals 1 outside fi, so this metric equals the Euclidean metric outside f2. Our 
goal is to reconstruct singularities of the initial data / given information about u on 
H xM. t . We therefore allow / to be any distribution in £'(SY). 

The solution to ([1]) can be expressed as a sum of two Fourier integral operators 
applied to /, as in |Duj : 



(3) u=(E + +£_)(/). 

E + and E- are of order —j with nonvanishing principal symbols, and have canon- 
ical relations Tl± which are, roughly speaking, the union of bicharacteristic strips. A 
bicharacteristic strip j(t) = (x(t),t,£(t),T) with initial condition (y,rj) is a smooth 
curve in T*(R™ x R) defined by 

d t (x(t),m) = ^(V e p,-X7 xP ) = (2c 2 (x)t,-V x (c 2 (x)M\ 2 ) 

(4) (z(0U(0)) = (y,v)- 

Since p is independent of t, r is constant and is determined by the condition p{pj) = 0. 
Thus given any initial condition (y, rf) there are two bicharacteristics with that initial 
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condition, which we will denote by 7 J ^(^,7^ Jt) where r > (resp. < 0) for 7+ 
(resp. 7~). We then have 

(5) n± = {{y,V,^t,^r) I {x,t,t,T)=~if yv) {t)} . 

Since c(x) = 1 outside fl, each bicharacteristic intersects H x R at most once for 
t > 0, and at most once for t < 0. 

Our goal is to reconstruct singularities of / from knowledge of the solution u on a 
portion d£l x R. However, we first consider a local problem in which we assume u is 
known on H x R. 

There exists T > such that the eikonal equation 

p(x, Vxfadtf) = \d t <jy\ 2 -c 2 (x)|Vx^| 2 = 

(6) 0(z,<U,t) = <*,0 

has two smooth solutions 0± on R" x [—T,T] t . We also have V x 4>± ^ on R™ x 
[— T, T], so since <j)± are homogeneous of degree 1 in £, we have for any K C R™ 
compact, 

(7) |V,0±(a;,t,OI>C|e| 

on K x [— T, T]. In fact, since c(x) = 1 outside f2, we can take C uniform on R n x 
[— T, T]. 

Dchne two Fourier integral operators S± of zcroth order by 

(S + f)(x',t) = jj e^'^-^a+ix', 0,t, V )f(y)dydr, 

(8) (S-f)(x',t) = ^^lle^'^-^a^x',0,t, V )f( y )d y d V . 

The amplitudes a± = a±(x,t,rj) G are constructed by geometrical optics so that 
the function 



(9) u(x',t):=(S + + S^)(f)(x',t) 

is the solution to (fTJ) restricted to H, modulo smoothing operators. This represen- 
tation is valid for \t\ < T. We have a+(x, 0, 77) = a_(x, 0, 77) = 1 and since they are 
constructed as solutions of transport equations, their principal symbols are nonvan- 
ishing. 

5+ and S- have formal adjoints acting on v G C^°(H x [— T, T]) given by 

(10) = III e^-^'^a.ix'AtMx'^dx' dtdv 

and they extend to continuous maps x [— T, T]) -> P(R n ). 

The canonical relations of S±, which we call H±, are defined when \t\ < T by 
taking a point (y, 77, a;', 0, t, £, r) G W± and projecting (a;',0, t, £,r) onto T*(x',t), a 
single fiber in T*(H x R) (see Figure 1). It is also possible to determine a point in 7i± 
given a point (x', t, r) G T*{H x [— T, T])\0 as follows: there is a unique number £„ 
such that p(x', 0, £ n , r) = and the bicharacteristic of p through (2', 0, i, £ n , r) 
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travels to the left (decreasing x n ) as t goes towards 0. Since f2 lies to the left of H, 
this is the only choice of such that the resulting bicharacteristic might lie in f2 
when t — 0. Continue this bicharacteristic through (a;',0,t, £ n ,r) back to {i = 0} 
and if it arrives in f2, project onto T*(f2) to obtain a point (y, rj). We can thus write 

(11) n± = {(y, n, x', t, C t) I (y, r), x', 0, t, £„, r) G H±} . 
From the local representation in ((9J we also have 

U± = {(V v </)±(x\0,t,T]),r),x\t,V( x , it )<f>±(x',0,t,r])) 

(12) | (x',i,r y )eIRr 1 x[-T,T] t xR; i } 
We note the following consequence of the preceding: 
Lemma 1.1. Suppose 

(13) V {x , tt) <t>±{x',Q,t,ri) = V( 9 ./ |t )^±(a/,0,t ) C) ■ 
TTien t] = (. 

Proof. Let £ n < be the unique number satisfying 

p(x',0, Va;'0±(a;',O,t, dt</>±(x',0,t, 17)) = . 

The existence and uniqueness of means the condition in (|13p determines a unique 
bicharacteristic strip 7 whose projection from T*(R™ xR t ) into M™ xt f passes through 
(af',0,t). Continue 7 back to {t = 0}. The projection of 7(0) from T*(R£ x R t ) into 
T*R™ is uniquely determined, so ([5]) implies that 77 = £. □ 

A point (y, 77) G T*(f2) will be called visible if there exists (x',t, £',r) such that 
either (a;', 0, t, £„, r) G 7^ or (x' ,0,t,£' ,£ n ,r) G 7 ( ~ r)) (with £„ determined as 
above). If \t\ < T, this is equivalent to having (y, 77, x' , t, r) G 7Y + or (y, 77, a;', f, r) G 
For example, if c = 1, (y, 77) is visible if and only if ?7„ 7^ 0. 

Consider the effect of applying S+ + SI to u. We must analyze the operators 
S+ 5+ and SIS-. Neither operator is well-defined on £'(£1), so cutoffs will have to be 
applied. Once this problem is fixed, it will turn out that the first operator reconstructs 
certain visible singularities of /, while the second operator is smoothing. 

In what follows, if / and g are distributions, we write / = g on U C T*(X)\0 if 
for any pseudodiffcrcntial operator A G L° whose essential support lies in U, we have 
A( f — g) G C°° . Two such distributions are called microlocally equivalent modulo C°° 
on U. 

Theorem 1.2. Let X G C$°(H x [— T, T]) and let V C H x [— T, T] &e an open 
set sttc/i t/ia£ x(x',t) = 1 /or aH (x',t) G Let C/± C T*(R") &e open cones in 
T*(R") smc/i t/ia£ i/ (y, 77) G L±, t/ie projection of the intersection of 7^ x ottj/i 
T*(H x [—T,T]) intersects H x [— T, T] in V 7 . TTien t/iere exisf pseudodifferential 
operators R± G L°(R") suc/i tfiai + )(S+ + SL)(/) = f onU+UU-. 

Proof. To analyze £+^5+ we follow the arguments of |GrSj| , chapter 10. We have, 
in the sense of distributions, 

(S* +X S+)(f)(z) = 4 ^ 2 ^ 2n y e 4(<2: ' c) -^ (a; '' ^' c)+<?i + (:l; '' '*'' )) - fe ' r ' )) 

(14) fl+ (i',0,t,C)a+(i'. 0, t, 77) X (x', t)f(y) dy d v dx' dt d( . 
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The integral is well-defined in the sense of distributions because x$+f £ £'(H x 
[-T, T]), and 5| :S'(Hx [-T, T]) -> £>'(R"). (We could approximate a+ by symbols 
in S* -00 to write S+X&+ as an absolutely convergent integral operator.) 
We denote the phase by 



(15) $ = (z, C) - <f>+{x', 0, «, C) + 4>+{x', 0, i, 7?) - (y, 77} ■ 

We break up the integral by inserting a cutoff -0 (^p) an d consider 



1 

4(2V) 2 



(S* +X S+)(f)(z) = — ^ / e'*« + ( I ',0,t,()a + (.T' I I ^) 



x V (^p) X(x\t)f{y)dyd v dx / dtd( 
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> a+(a;',0,t,C)a+(a;',0,i,77) 



4(2tt) 2 ™ 

(16) x ^1 _ ^ r liz^l n x (a/, t )f(y) dy d V dx> dt d( . 

By Lemma 11.11 $ has no stationary points in the support of the amplitude of the 
second integral, so we only need to consider the first. Write r\ = Xto where A = \rj\ 
and make the change of variables £ = Act. Then we are left with 

(17) -^-^ J e i ^a + (x',0,t,Xa)a + (x',0,t,r 1 )^(\a-u\) X (x',t)f(y)dydr,dx'dtda 
where 



(18) 1= (z,a) - <f>+(x',0,t,a) + <t>+(x',0,t,uj) - (j/.w) . 

We use the method of stationary phase on the (x' , t, a) -integrals. (The integrand 
is compactly supported in these variables.) By Lemma ll.ll we have the critical point 
ct = lj, z = V cr 0+(a; / , 0, t, cr)\o-=u- In a neighborhood of any point of H.+ , H.+ is 
parametrized by the initial conditions (y, rf) of bicharacteristics, so TL + is a smooth 
manifold of dimension 2n. This implies that the critical point is non-degenerate (see 
[Tre] . Proposition 8.1.2). 

(S* +X S + )(f)(z) = ^-JJ e ^-y^b(z, V )dyd V 

(19) + Kf(z), 

where K G L~°°. b 6 Sf and the principal symbol is 

(20) b (z,r)) = \a + {x' ,Q,t,rj)\ 2 x{x' ,t) \v v <p+(x> ,o,t, v )=z ■ 
In light of (fT2|) . we can write 



(21) b {z,T]) = |a + (7roH + (z,r;),r / )| 2 x ( x o tt o H+)(z, ry) 

where 7r : T*(H x [—5, 5]) — > i? x [— T, T] is the natural projection. 

60 7^ on U+, so there exists a pseudodifferential operator G L°(U+) such that 
(Rf + S +X S+)f = fon U+. 



By a similar procedure, we can also construct R'_ € L°(U-) such that (R'_S~xS~)f 
/onV. 

We next examine the operator S+xS-- We have 

(S+xS-)(f)(z) = ~ J e i (( z >O~'l>+{x',0.tX)+<P-(x',0,t,ri)-{y,r))) 

(22) a + (x',0,t,C)a-(x', 0, t, rj)x(x', t)f(y) dy dt] dx 1 dt dC, . 
so the Schwartz kernel of SI xS- is 

(23) K(z, v) = J e^^-^'^^+^-^'^^-^^Aix', t, rj, Qdf] dx' dt d( , 
where 



(24) A(x',t,ri,0 = ^r^a + (x',0,t,Oa-(x',0,t,r))x(x',t) 
is of class Si . 

The derivative in t of the phase could vanish only if dt4>+ — dt<f>-. But <j) + and 
0_ are distinct solutions of the eikonal equation, so their t derivatives have differing 
sign. In fact, 

d t (j)+{x,t,ri) = c(x)\V x (/) + (x,t,f])\ 

(25) dt4>-{x,t,rf) = -c(x)\V x 4>-(x,t,T))\. 
Define a differential operator 

(26) L = i(d t <t>-(x',0,t,r,) - d t <l> + (x',0,t,O) dt • 
Then L fixes the exponential portion of the integral: 

(27) L (^ e i ({ z X)~<P+(x',tX)+4>-(x',t,r,)-(y,ri))^ _ e i({z,C) -<t>+ (x',t X)+<P- ,t,v) - (v,V)) _ 

From Q we have 

(28) ,,, r 1 ^ 7-r^ < ( 



\dt4>- {x 1 , 0, t, rf) - d t (j> + (x', 0, t, C) | " |??| + |CI ' 
Since dtt<fi±(x, t, £) is also homogeneous in £ of degree 1, we have 



(29) 



d t ^_(a:',O,t,» 7 )-d t 0+(a:',O,t,C) 



C 

< 



M + ICI 



After repeated integration by parts using the operator L*, we get 



K(z,y) = I e <«*.0-*+(* , ,o,t,0+*-(* , .o,t,u)-<»,n» {{L t ) k A{x',t,rj,0)dr 1 dx'dtdC > 
(30) < / e 4(<z ' C> -^+ (:E '' ^' C)+<?i - (:l; '' ^' T ' ) ~ ( ^ ? ' )) A_ fc (a;',i,77,C)d?7dx' AdC 
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where A_ k G S 1 q . This can be achieved for any k > 0. Hence K G C°°(R 2n ), and 
S+\S- € A similar argument also shows 5* \S+ G 

We now combine i?^ and to form a full parametrix for x(S+ + S-). Let 
^±(2/7 ^) £ C°°(E/+ U [/_) be homogeneous functions of degree forming a partition of 
unity on U±. Let 

(31) R±{x,D) = 9±(x,D)R' ± {x,D) . 

Then by construction, on £/+ (1 [/_ we have 



( J R + S; + J R_^)(S+ + £_)(/) 

(32) 



(0 + ^S;S + )(/) + 9-R'_S*_ SL)(/) 

e+f+e-f 
/■ 

□ 



If the metric has no conjugate points, then we can take T = 00. As a result, 
for any visible singularity (y, 77) we can choose V such that (y, rf) lies in U+ U [/_. 

2. Reconstruction of the Wavefront Set Given Measurements on the 

Boundary 

We now investigate the situation when the measurement surface L is not a hyper- 
plane, but instead is a relatively open subset of d£l. For simplicity we assume that f2 
is strictly convex, so that any bicharacteristic along which a singularity might travel 
intersects dQ. at most once for t > and also at most once for t < 0. Let T be an open 
subset of dfl compactly contained in F. We will attempt to reconstruct singularities 
that reach F. 

Let W\ , . . . Wk be a partition of T such that we have boundary normal coordinates 
near each Wj. Then, as before we define Uj.± C T*(K") be open cones in T*(R n ) 
such that if (y, 77) G f/j,±, the projection of the intersection of 7^ * with T*(dCl x ffi.) 
intersects 3!lxM in Wj x [-T.T]. Let X j € Co°(T*(dfi x E))'be cutoff functions 
such that supp(Yj) C Wj X K and JZjXj = 1 on T X [— T, T]. Our measurements 
correspond to knowing w|f x j_ tt j. Let 

(33) E/ = U^,±- 

j,± 

We call U the visibility set, and attempt to reconstruct the singularities of / on U. 

Consider a single patch Wj and let (w',w n ) be boundary normal coordinates near 
Wj such that fl corresponds to w n < 0. Writing Xj{S+ + S-) ul those coordinates 
we obtain data in the same form as in ([5]): 

Uj (w',t) = + £*_)(/)(«/,<) 

(34) + — fj e^-Vtri-lwVaj,-^ 
where 
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ip jt ±(w',t,ri) = <j>±(x,t,ri) 

(35) a ji± (w',t,ri) = a±{x,t,rf) 

on Wj x [-T,T]. 

We apply the operator 

(Sl + + Sl_)(.)(z) = HI e^-^'^-a^'^V^dw* dt d V 

(36) + HI e^-^-^'^d^'^v)^) dw' dt dr, 

to Uj. As above, we obtain pseudodifferential operators of order 0, modulo smoothing 
operators, which are elliptic on Uj. + and Uj-, respectively. Hence there exist Rj,± G 
L\ elliptic such that 

(37) (Rj, + S*, + + Rj,Sl_) Xj (S j!+ + Sj,-)f = f 
on Uj. + U Uj-. 

Let 6j = 0j(y,r}) G C°°{U) form a partition of unity on U with each 0j supported 
on Uj t + U Uj t -. Then we have 

(38) Yl 9 ^ D y)( R i,+ s l+ + R j,- s l-)Xj(Sj,+ + Sj,-)f = f 

3 

on U. 

In conclusion, we discuss a simple condition which guarantees that we can reson- 
struct the entire wavefront set of /. Suppose that for each (y, n) G T*(Q), cither 7 f . 

or 77 , reaches T in a time less than T. (Equivalently, we may say that either of the 

unit speed geodesies eminating from (y, rf) or (y, —rj) exits fl through T in a time less 
than T.) Then every singularity of / is visible, i.e. T*(Q) CU. As a result, 

(39) X 6 i(v> D y ){R 3 , + Sl + + Rj,-Sl_) X j 

3 

reconstructs / from u\f, x ^ TT y We summarize this discussion in the following 
theorem: 

Theorem 2.1. Let f e £'(Q) and Zet U be the visible set. Then with Xj, Sj,±, Rj,± 
as above, 



(40) 
on U. 



£ 0j(y, D y ){R h+ Sl + + J? j ,_5*_) Xj (% + + Sj,-)f = f 
3 
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Remark 2.2. This work is part of the author's dissertation at UCLA. Recently a 
paper [StUhlj by P. Stefanov and G. Uhlmann was announced in which similar results 
were obtained. Our proofs differ from those in [StUhlj . 
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